An Uncertainty Principle for Ensembles of Oscillators Driven by Common Noise 
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We discuss control techniques for noisy self-sustained oscillators with a focus on reliability, stabil- 
ity of the response to noisy driving, and oscillation coherence understood in the sense of constancy 
of oscillation frequency. For any kind of linear feedback control — single and multiple delay feedback, 
linear frequency filter, etc. — the phase diffusion constant, quantifying coherence, and the Lyapunov 
exponent, quantifying reliability, can be efficiently controlled but their ratio remains constant. Thus, 
an "uncertainty principle" can be formulated: the loss of reliability occurs when coherence is en- 
hanced and, vice versa, coherence is weakened when reliability is enhanced. Treatment of this 
principle for ensembles of oscillators synchronized by common noise or global coupling reveals a 
substantial difference between the cases of slightly non-identical oscillators and identical ones with 
intrinsic noise. 

PACS numbers: 05.40.-a, 02.50.Ey, 05.45.Xt 
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Collective phenomena in ensembles of dynamical 
systems can manifest complex behavior and self- 
organization, which are not possible for unitary systems 
with arbitrary level of complexity. For biological sys- 
tems, they extend from the simplest collective dynamics 
of bacteria [l[ or higher organisms 0] to the activity of 
neural tissue [1, 0] ■ Similarly, in technology, essentially 
collective phenomena vary from plain synchronization ef- 
fects [1, @| to the phenomena laying in the basis of oper- 
ation of artificial neural networks [3, [3| . Although these 
phenomena are inherent to large ensembles, their features 
are determined by characteristics which make physical 
sense also for unitary systems @. Indeed, synchroniza- 
tion (and clustering) in ensembles and their susceptibility 
to control are influenced by their individual robustness 
properties. "Reliability" (or "consistency" [9j), i.e. 
the stabilityof the system response to noisy driving, and 
coherence [13], i-e. the constancy of the oscillation 
frequency in time, are the principal ones. 

Mathematically, the natural quantifier of reliability is 
the Lyapunov exponent (LE, A) measuring the expo- 
nential decay rate of perturbations of the system re- 
sponse [1, [ll|, 12 1 . Oscillators are more reliable for a 
larger negative LE. Coherence can be quantified by the 
diffusion constant (DC, D) of the oscillation phase (p{t), 
which grows non-uniformly in time due to either noise or 
chaotic dynamics; {{(p{t) — {if)t)'^) oc Dt. These quan- 
tifiers are immediately relevant for the synchronization 
phenomenon in ensembles of oscillators subject to com- 
mon noise [l3, [3) or global coupling [1, [l^ 13] • 

In this Letter, we report how the reliability and co- 
herence of a noisy limit-cycle oscillator can be efficiently 
controlled by a general linear feedback. However, the 
ratio of LE and DC is shown to be independent of the 
noise strength and feedback parameters. The persistence 
of this ratio can be formulated as an uncertainty princi- 
ple. Finally, the important implications of this principle 
for controlling collective dynamics of ensembles are re- 



vealed, and conclusions are drawn. 

Let us consider an A^-dimensional limit-cycle oscillator 
subject to small linear delayed feedback and weak noise: 



X, =F,(x)+az,(t) + B,(x)o^(t), 



(1) 



where i = 1, 2, . 
back term 

Z^{t) 



,N, ais the feedback strength, the feed- 
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Gij{t) is the Green's function, different Gij{t) can fea- 
ture single or multiple time delayed linear feedback (par- 
ticular forms of this function are specified in the text 
below) [l3,[l3i linear frequency filter etc.; "o" in- 
dicates the Stratonovich form of equation, £^{t) is white 
Gaussian noise: = 0, {£.{t)£,{t')) = 25{t - t'). 

For weak noise and feedback the dynamics can be de- 
scribed within the framework of the phase description up 
to the leading order of accuracy [3, HO] : 
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b = no + a I Y,Y.G,jih)H,,{if{t-h),ip{t))dh 
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+£/(<^(i))°CW, (2) 
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where i7o is the natural frequency of the oscillator, / 
is a 27r-periodic function featuring the sensitivity of the 
phase to noise, e is the noise amplitude, Lp) is the 

increase of the phase growth rate created by the feedback 
term Xj{'ip) acting on the variable 2:^(1^9). The phase of 
the noise-free oscillator grows uniformly and its shifts nei- 
ther decay nor grow in time. Noise creates irregularity of 
the phase growth rate, measured by the phase diffusion 
constant D (DC): {{ip{t) - {ip{t))f) cx Dt. Addition- 
ally, noise results in convergence of trajectories, phase 
shifts decay, and the exponential rate of this decay is 
measured by the Lyapunov exponent A (LE). Laborious 
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FIG. 1: (Color online) Dependencies of the Lyapunov exponent A and the diffusion constant D (upper and lower graphs in 
the upper row of plots, respectively) on the feedback/filter parameters for the van der Pol oscillator with fi = 0.7 subject to 
Gaussian white noise of strength e = 0.05. (a) "Simple" delayed feedback [Eq. (|8])]. (b) Multiple time delayed feedback [Eq. ((9])]. 
(c) Linear frequency filter [Eq. (|1U|) ] with a = 0.1. The oscillation period of the control-free noiseless system Tb « 27r/0.97 ; the 
feedback strength (a = 0.03 and a = —0.03) is specified in plots. Solid black lines present analytical dependencies [Eqs. ©-((S])] 
with (/^) f» 0.5464 and ((/')^) ~ 0.5775. The ratio between the Lyapunov exponent and diffusion constant is constant and 
obeys Eq. ((6]) up to the calculation accuracy (see the lower row of plots). 



but straightforward evaluation (non-general versions of 
which can be found in EEHEil]) yields 
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D = 
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(4) 
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(5) 



where (...)^ = (27r) ^ J^"" ...dip, htj{->p) := {H^j{ip + 
^p,(p))^, and the prime denotes derivative. Notice, for 
any kind of linear feedback the ratio 
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is independent of parameters of the feedback. 

The modified formulation of the result ([B]) was revealed 
also for a certain class of systems where the phase can- 
not be well defined. Specifically, in systems below a Hopf 



bifurcation, noise can induce oscillatory motion [22*1 for 
which the phase is not well defined. For slightly per- 
turbed periodic oscillations the autocorrelation function 
Cjj{s) := {xj{t)xj{t + s)) decays exponentially, Cjj{s) = 
Cj{s) exp(^D\s\), where function Cj{s) oscillates but nei- 
ther grows nor decays asymptotically. One can introduce 
an alternative quantifier of the irregularity of oscillations, 
the correlation time icorr := C~^(0) \Cjj{s)\ds and 
find tcorr oc l/D for perturbed periodic oscillations. For 
noise-induced oscillations, icorr is an appropriate quanti- 
fier, while D cannot be evaluated. In Refs. for the 
van der Pol oscillator slightly below a Hopf bifurcation, 
the product icorrlAI was reported to be independent of 
the feedback with time-delay (s), although the feedback 
significantly changed icon- and A. This is equivalent to 
the result (©. 

The validity of our findings can be underpinned 

with the results of numerical simulation for noisy van der 
Pol oscillator 



y, 



y ^ ^{1 - Ax'^)y - X + az{t) + e ^{t) . (7) 



Here describes closeness to the Hopf bifurcation point, 
the phase (p — ~ arctan(y/a;) (for ^ ^ 1 the limit cycle 
is: X — cosip, y — ~ siiup). We consider 3 possible cases: 
(1) a "simple" delayed feedback with delay time r; 



z(i-)(t) = 2(y(t-T)-y(i)), 



(8) 
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FIG. 2: (Color online) Ensemble of 100 van der Pol oscil- 
lators ((7| with nonidentical frequencies subject to the sim- 
ple delayed feedback The distribution of frequencies is 
the gaussian one centered at 1 with standard deviation 0.001, 
feedback strength a — 0.03, noise amplitude e = 0.1. (a) The 
gray line plots trajectory of one oscillator, the symbols plot 
snapshots of the ensemble for two specified values of the delay 
time (check these values in Fig. [T^). (b) The distribution of 
the phase deviations from the value ipo corresponding to the 
instantaneous ensemble-mean state, (c) The deviation of the 
phase from its mean growth for an oscillator (the offset of the 
vertical axis is arbitrary). 



which yields 

G(i")(t) 





2{S{t - t) ~ 6(t ^ 0)) 



(here we exphcitly indicate that J^°° 5{t — 0)dt = 1). 
(2) a multiple time delayed feedback; 
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{t)^2j2 {y{t - + 1)^) - y{t - nr)) , (9) 



where \R\ < 1, which yields 



G("^)(t) 




2j^R''{d{t-{n + l)T)-d{t-nT 
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(3) a linear frequency filter; 

z'^^^^t) ^2u{t), u + au + ujlu = ax{t), (10) 
which yields 
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FIG. 3: (Color online) Ensemble of 100 van der Pol oscillators 
([7]) subject to the simple delayed feedback (|8]), common noise 
£ = 0.1 and intrinsic noise eint = 0.005. For details see caption 

to Fig. [a 



For /i ^ 1, the van der Pol oscillator possesses a cir- 
cular limit-cycle and h2i{ip) = — (l/2)cos(^, /i22(<y3) — 
(1/2) sin f{(p) = —sinLp. The analytical results de- 
rived with the above formulae match well the results of 
lol . 17 . isl . as can be well seen in 



numerical simulation 
Fig. [T] Remarkably, the ratio between LE and DC is con- 
stant with a good accuracy even when they vary by one 
order of magnitude, in agreement with Eq. ([6|). 

Constancy of the ratio ^ can be formulated as a kind 
of ^^uncertainty principle^^ for the linear feedback control 
techniques: 

The reliability of a noisy oscillator can be significantly 
enhanced (by means of a weak linear feedback) but at the 
price of the loss of its coherence and, vice versa, the co- 
herence can be significantly enhanced but with the loss of 
reliability. 

For ensembles of uncoupled oscillators driven by com- 
mon noise we know, that the imperfectness of iden- 
tity of oscillators or intrinsic noise leads to imperfect- 
ness of synchronization. The characteristic spreading of 
states (X 1/^/^^ [14,] • That is, higher reliability leads 
to stronger synchrony of oscillations. Thinking of em- 
ployment of control techniques we should distinguish two 
"pure" situations: (i) slightly non-identical oscillators 
driven by identical noise and (ii) identical oscillators with 
small intrinsic noise, individual for each oscillator. 

For an ensemble of non-identical oscillators, the dis- 
persing of phases is owed to the mismatch of natural fre- 
quencies. This mismatch is nearly unaffected by small 
feedback, while a synchronizing action measured by LE 
is influenced by the feedback. Hence, the dispersion of 
states is perceptive to the feedback control. For instance, 
one can see in Fig. [2] that the ensemble is well syn- 
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chronized (Fig. [2^,b) when the phase diffusion is large 
(Fig. [2}:), and, on the contrary, for a smaU diffusion (en- 
hanced coherence) we observe a much poorer synchrony. 

For an ensemble of identical oscillators with intrin- 
sic noise, the situation is different. The dispersion is 
owed to the mutual diffusion of oscillator phases cre- 
ated by intrinsic noise. Similarly to we can find 
Ay) cx VMnt/V— A, where Dint characterizes the diffu- 
sion due to intrinsic noise. The diffusion owing to in- 
trinsic noise is expected to be subjected to the effect of 
the feedback in the same manner as the total diffusion. 
Hence, we expect the feedback to not influence the phase 

const. Indeed, in Fig. [3] 



dispersion, Aip cx V-Dint/ v— A : 
one can see that the distribution of phase deviations is 
tolerant to feedback control. In this case, one can employ 
the feedback control techniques to control the coherence 
without the loss of synchrony. 

Summarizing, we have discovered the universal uncer- 
tainty principle, which is valid for general noisy limit- 
cycle oscillators subject to a general linear feedback con- 
trol and proved it both analytically and numerically. 
Mathematically, this uncertainty principle takes the form 
of Eq. (O; the ratio of the Lyapunov exponent (A), mea- 
suring reliability (or the response stability, [1, Q), and 
the diffusion constant (D) , measuring coherence ( or the 
constancy of the instantaneous frequency in time), is in- 
dependent of the noise strength and feedback parame- 
ters. That is the reliability of a noisy oscillator can be 
significantly enhanced by means of a relatively weak lin- 
ear feedback, but at the price of the loss of its coher- 
ence, and vice versa, the coherence can be significantly 
enhanced but with the loss of reliability. The principle 
has an implication to practical issues of the control of 
synchronization in non-ideal ensembles of oscillators in 
nature and technology [J, Q, 0, 0, [3, For ensem- 

bles of weakly non-identical oscillators driven by common 
noise (cf Fig. [2|) the enhancement of synchrony, achieved 
due to the reliability increase, results in a poorer coher- 
ence of each individual oscillator. Unsimilarly, for ensem- 
bles of identical oscillators with intrinsic noise (cf Fig. [3]) 
the synchrony is not influenced by the change of relia- 
bility/coherence; therefore, coherence can be controlled 
without stray effects on the synchronization. 

Notice, that ensembles of globally coupled oscillators 
in an asynchronous state immediately correspond to the 
case of uncoupled oscillators receiving common driving; 
this driving is the ensemble-mean value forcing an os- 
cillator via the coupling term [15;]. As long as collec- 
tive modes vanish, the ensemble-mean value fluctuating 
about zero may be well considered as a stochastic pro- 
cess and one can speak of synchronization by common 
noise. (When the collective mode appears, it can be reg- 
ular in time; synchronization by a regular signal drasti- 

■) 



lators as well [16 1 
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cally differs from that by a nonperiodic one, e.g. 
Therefore, our findings are related to the fundamentals 
of synchronization in ensembles of globally coupled oscil- 
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